The Quantum Torus Chain 
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We introduce a new set of one dimensional quantum lattice models which we refer to as The quan- 
tum torus chain. These models have discrete global symmetry, and projective on-site representations. 
They possess an integer- valued parameter which controls the presence or absence of frustration. De- 
pending on whether this parameter is even or odd these models either exhibit gapped symmetry 
breaking phases with isolated critical points, or gapped symmetry breaking phases separated by 
gapless phases. We discuss the property of these phases and phase transitions for two special values 
of the parameter and point out many open problems. 



I. INTRODUCTION 

Over the years many interesting lattice models have 
been introduced to capture the essence of important 
physical concepts, and make them open for more quan- 
titative studies. In modern language quantum lattice 
models can capture short-range entangled states such as 
symmetry breaking or symmetry protected topological 
states 1 -. It also can capture long-range entangled states 
such as quantum critical and the "intrinsic topological 
ordered" states. 

The importance of symmetry is well known in classi- 
cal and quantum physics. Recently it is realized that 
a strong tie exists between quantum entanglement and 
symmetry^ For example, in one space dimension Chen 
et a£ showed that short-range entangled, fully symmet- 
ric, quantum states are classified by the projective repre- 
sentation of the internal symmetry group. In addition, 
if the site representation of the symmetry group is pro- 
jective, short-range entanglement is impossible without 
symmetry breaking. Because in one dimension a gapped 
system is necessarily short-range entangled, this implies 
projective on-site representation can not have an energy 
gap without breaking some symmetry. This result is a 
generalization of Haldane's work on SO (3) spin chair* 2 -. 
There half integer spin chains need to break a symmetry 
to open an energy gap, while integer spin chains do not. 

In this paper we introduce a new family of one di- 
mensional quantum lattice models whose global symme- 
try groups are discrete and the on-site representation is 
projective. Indeed, these models are either gapless or 
exhibit spontaneously symmetry breaking. In addition, 
they possess an integer- valued parameter which controls 
the presence or absence of frustration. Depending on 
whether this parameter is even or odd, qualitatively dif- 
ferent phase diagrams are observed. By applying the 



density matrix renormalization group (DMRG^ and 
Matrix product state (MPS)£ methods, we find for even 
parameter these models are generically gapped and show 
spontaneous symmetry breaking. Fine tuning is required 
to close the energy gap. On the other hand for odd 
parameter, both gapped symmetry breaking phase and 
gapless symmetric phase are generic. The results are all 
consistent with the conclusions drawn by Chen et. a&. 

The results reported in this paper show that the new 
class of models that we introduced has a rich set of inter- 
esting phases. We hope that the study of this new class 
of models will enhance our understanding of symmetry 
protected topological phases in one dimensional systems, 
and of topological phases in a more general setting. 



II. THE QUANTUM TORUS CHAIN 

The family of models, which we refer to as the quantum 
torus chain, is motivated from the following view of the 
quantum spin chain. 

The dimension D of the single site Hilbert space of a 
SU(2) spin-5 1 chain is 25+1. The basis states can be 
viewed as the one-particle states of a unit-charge parti- 
cle running on a sphere enclosing a magnetic monopolc 
which produces 25* magnetic flux quanta through the 
sphere^. The nonzero total Gaussian curvature of the 
sphere leads to a term in the Hamiltonian which has 
the same form as the term describing the flux. Thus, 
the curvature effectively increases the flux, and there are 
D = 2S + 1 rather than 2S one-particle states (see also 7 -). 
The spin operators S x,v,z are the generators of the mag- 
netic translations. 

In the XYZ model the nearest neighboring spin opcr- 
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(a) 




FIG. 1: (a) Each site of a quantum spin- 5 chain can be viewed 
as the lowest Landau level of a unit charged particle moving 
on a sphere under the magnetic field of a Dirac monopole. 
The total magnetic flux produced by the monopole is 2S. (b) 
The Torus chain is a modification of the spin chain where each 
sphere is replaced by a torus. 



ators are coupled as 
Hxyz = / ](JxSfS? +1 + J y SfSf +1 + J z SfSf +1 ). (1) 

i 

This is depicted in Fig. HJa). 

In this work we consider the "torus chain" (Fig.QJb)), 
where the spheres of Fig. [lja) are replaced by the 
torus. Because the total Gauss curvature of a torus 
is zero, the lowest Landau level degeneracy is equal to 
the number of flux quanta through each torus^. Un- 
like the sphere, the magnetic translation group of a 
torus, the Heisenberg-Weyl group, is discrete rather than 
continuous^. The commutation relation between the op- 
erators U x , U y which generate the magnetic translation 
in two orthogonal directions is 

U X U V = e l2 * /m U v U x , (2) 

where (U x ) m — (U y ) m = I. The m in the preceding 
equation is the number of flux quanta through each torus, 
and is the integer valued parameter we referred to earlier. 
In the following we shall encounter two different bases, 
\q) and \q), for the site-Hilbert space. Both q and q are 
defined modulo m. In the |g) basis, U x and U v are given 
by 

m—l 

U x = ^ q/m \q)(q\, 

rn — 1 

tf*=X> + l)(g|. (3) 
The dual basis \q) is the Fourier transform of |g), i.e. 

-. m— 1 

\q) = ^ V e 2 ^ m \q) . 

In this basis 

m—l 

u x = j2 i?+i><?i, 

9=0 
m—l 

Uy=J2^ i9/m \Q)(q\- (4) 

q=0 



Analogously to the XYZ model we couple the gener- 
ators of the magnetic translation group of the torus to 
form the following Hamiltonian: 

H = ]T(cos0 U x U x l x + sin0 UfUfl, + h.c.). (5) 

i 

Here, the parameter 6 is introduced to control the relative 
strength (and sign) between two terms, J^i UfU'^ + h.c. 
and J^i + h.c, in the Hamiltonian. In the follow- 

ing we assume the number of sites L to be an integer 
multiple of m. With this restriction the two unitary op- 
erators U x = n£=i Uf and U v = Y\ i=1 Uf commute with 
each other and with the Hamiltonian ([5]). The Hamilto- 
nian is one of the simplest models which preserve the 
full group of magnetic translations (the group generated 
by U x and W\ i.e. an analogue of an SU(2) symmetric 
spin chain. It is possible to consider a further generaliza- 
tion of the model, with the same symmetry, by adding 
different complex phase factors for the two terms. How- 
ever for simplicity we shall restrict ourselves to the case 
where the coupling ratio is real. We also note the simi- 
larity of our model to the m-state Potts model (although 
the symmetries are a bit different). For more descriptions 
about related/similar models, see, e.g., Refs. I9j-fl2l. 

We will denote the conserved quantum numbers asso- 
ciated with U x and W> by ( e 27ri( 3/ m , e -2niQ/m^ or morc 
compactly as (Q, Q). In terms of qi and qt we have 

Q = qi mod m , Q = qi mod m . 

i i 

In addition to IA X and U y there are other symmetry op- 
erators that leave the Hamiltonian ([S]) invariant. First of 
all, there is the g-invcrsion, 1Z = Yii=i Ri > where 

RMi) = I - Qi) mod m, Ri\qi) = |— g,-) mod m. 

Secondly, there are the following anti-unitary operations, 
K and K = RK, which cause complex conjugation in the 
g and g-basis respectively. We note that K 2 = K 2 = I. 
The operators U x , U v , 1Z, and K generate the internal 
symmetry group G. As one can readily check the group 
multiplications of G are identical to those of the symme- 
try group of a rectangular, periodic lattice on a torus, 
with U x and U v as the lattice translations, R as rotation 
of angle it, and K and K as the two reflections. We note, 
as usual, the translation generators of the space group 
commute. This is because we have restricted L to inte- 
ger multiple of m. It is also interesting to note, unlike 
the usual cases, some of the space group elements are 
represented anti-unitarily here. The internal symmetry 
group G extended with the translations and inversion of 
the one-dimensional (ID) chain defines the full symmetry 
group of the Hamiltonian. In addition the Hamiltonian 
([5]) has a "duality" symmetry upon «-> tt/2 — 8 (mod 2ir) 
and U x <-> U y . 

The irreducible representations of the group of internal 
symmetries determine the (minimal) degeneracies of the 
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energy levels. For example for m = 3, which we shall 
spend great length discussing, the multiplets take the 
following form in terms of the (Q,Q) 

(Q,Q) -(0,0) 

(Q,Q) -(1,0); (2,0) 

(Q,Q) =(0,1); (0,2) 

(Q,Q) =(1,1); (1,2); (2,1); (2, 2) 

At the "self-dual" point tan 9 = 1 the two terms in the 
Hamiltonian ([5]) have the same coefficient. There the 
symmetry group is larger, which causes the two doublets 
in the above table to become degenerate, and form a 
quadruplet. 



(a) (b) 




FIG. 2: The phase diagrams of Eq. © for (a) m = 2 and (b) 
m = 3. The black and red dots mark first order transitions 
and critical points, respectively. At points marked by crosses 
the model is classical hence exactly solvable. 



III. SOME PROPERTIES FOR GENERAL m 

Having introduced the model and its symmetries, we 
start our analysis by first considering the dependence of 
the model on the parameter to. It turns out that even 
and odd to have qualitatively different phase diagrams. 

For later discussions it is useful to consider the sub- 
group Q generated by the subset of generators that are 
unitary: {I,U X , U y ,lZ}. This group is isomorphic to the 
symmetry group of an to x to oblique lattice on a torus. 
It is important to note, however, for each site of the ID 
chain this symmetry group is represented by a projective 
representation. The non-trivial U{\) phases are due to 
the presence of the factor e i2lT / m in U X U V = e i2 V™ [/!/[/*. 
Applying the results of Chen et al we therefore conclude 
that Eq. (O should cither exhibit a gapped spectrum with 
a symmetry breaking ground state, or a gapless spectrum 
with a symmetric ground state. 

To see how this is realized, we first consider the sim- 
plest case, with m = 2. In this case we can rewrite Eq. ([5]) 
in terms of a spin 1/2 Hamiltonian. This is achieved by 
identifying q = 0, 1 with spin up and down, and express 
U x ' y in terms of the Pauli matrices U x — > a z ,U y -> a x . 
Under this identification Eq. ([5]) becomes 

H -> 2 cos 9J2 + 2 sin 9 £ of af +1 . 

i i 

This is the anisotropic XZ model, which is gapped for 
all 9 except 9 = ±7r/4 and ±57r/4 where it is quantum 
critical. The phase diagram, which is symmetric under 
9 — > 7r/2 — 9 due to duality, is shown in Fig. [U[a). In the 
gapped phases the U x jU y and/or the translation sym- 
metry is spontaneously broken leading to two-fold ground 
state degeneracy. At the points 9 = ±7r/4 and 9 = ±57r/4 
the ground state respects all symmetries but the energy 
spectrum is gapless. 

To get a qualitative understanding of how this phase 
diagram generalizes to other values of to, it is instruc- 
tive to consider four special points. First, for 9 = tt 



the Hamiltonian becomes classical because there are no 
non-commuting operators. In that case the ground state 
is given by Y\f=i W) an d is to- fold degenerate for any 
given to. Clearly the U y symmetry is spontaneously 
broken and there is an energy gap. The gap persists 
for small 9 — n, which shows that this is an extended 
gapped phase. For to = 2 it extends to the full interval 
37r/4 < 9 < 57r/4, with an equivalent, dual phase in the 
interval 57r/4 < 9 < 7-k/A. 

Provided to is even the situation is similar at the point 
9 = 0. In that case the ground state is also to- fold degen- 
erate, now with a ground state of the form Yli=i \ (~^Y l)- 
Also around the point 9 = there is a gapped phase, 
with a ground state that, in the L — > oo limit, spon- 
taneously breaks both IA V and translational symmetry. 
For to = 2 this phase also extends to the full interval 
— 7r/4 < 9 < 7r/4, with the equivalent, dual phase in the 
interval 7r/4 < 9 < 37r/4 

However, if to is odd the situation is quite different. 
For 9 = 0, all states with 

q% - q%+i = (to ± l)/2 (mod to) 

yield the same energy. As the result, the ground state 
is extensively degenerate. For example, in an open chain 
the degeneracy is 2 i_1 to. When 9 deviates from zero 
these degenerate configurations are mixed by the sin# 
term, and hence the degeneracy is lifted. This is very 
similar to the effects of quantum fluctuations in frus- 
trated magnets. Due to the duality of the Hamiltonian, 
the same discussions hold for 9 = n/2, except the roles 
of q and q are interchanged. 

IV. THE CASE m = 3 

In the previous section, we used symmetry considera- 
tions to gain knowledge about the behavior of the model 
at some special points, and some general properties of 
the phase diagram. In this section, we will focus on the 
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case m = 3. We start by giving a quick overview of 
the phase diagram in the next subsection, followed by a 
more detailed description of the various phases and phase 
transitions in the subsequent subsections. 



A. Overview of the m = 3 phase diagram 

We hrst show the exact diagonalization results of a 
small system (namely, L = 12 sites). In Fig. [31 the 

Ground state energy per site, L=12 
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FIG. 3: Upper panel: ground state energy per site of the 
m — 3 model for L = 12 sites, as a function of the angle 6. 
Lower panel: first derivative of the ground state energy per 
site of the m = 3 model for L — 12 sites, as a function of the 
angle 6. 



ground state energy per site is shown in the upper panel, 
while the first derivative is given in the lower panel. 

The figures give a clear indication of a first order tran- 
sition, due to a level crossing, at 6 = 0. This is sig- 
naled by a sharp kink in the ground state energy. As 
already discussed, there is at this point also an extensive 
ground state degeneracy, which grows exponentially with 
the system size. Due to the 'duality symmetry', there is 
an analogous first order transition at 9 = 7r/2. 

This first order transition connects two gapless phases. 
In fact, we will later show that throughout the whole re- 
gion < 9 < tt/2 the system is critical and characterized 
by a central charge c = 2. This region contains the special 
point 9 = 7r/4, where the Hamiltonian is self dual. More- 



over we will show in the next subsection that the critical 
behavior for 9 < can be related to that for 6 > (and 
similarly the behaviour for 9 > tt/2 is related to that of 
9 < tt/2) . The critical region for 9 > tt/2 gives way to 
a gapped phase at around 9 ~ 0.67T, which is signaled by 
a sudden drop in the first derivative of the ground state 
energy, see Fig. [3J A much more detailed study of this 
phase transition will be given in Subscc. IIVDI 

Beyond 9 w 0.67T, and in fact in the whole region 
0.67T < 9 < 5tt/A, the system is gapped, with a three- 
fold degenerate ground state in the thermodynamic limit. 
In this gapped phase the IA V symmetry is spontaneously 
broken. The property of this phase can be understood by 
considering the special point 9 — tt, where the Hamilto- 
nian has no non-commuting operators and can trivially 
be solved. 

The gapped phase in the region 0.67T < 6 < 5tt/4 has 
its dual analog in the region 57r/4 < 9 < 1.9n, where 
the 1A X symmetry is spontaneously broken, giving rise to 
a three-fold degenerate ground state. In addition to the 
above the two panels in Fig. [3] also give hints of a phase 
transition between the two gapped phases at 9 = 5tt/4. 
Subsec. IIVFI is devoted to more detailed discussion on 
this phase transition. 

B. The behaviour near 6 = 

Next, let us focus on 9 w 0, and for simplicity let us 
concentrate on m = 3. The arguments given below can 
be applied for arbitrary odd m. As discussed above the 
ground states of the cos 9 term of the Hamiltonian con- 
sists of |{<Zj}) configurations where no nearest neighbors 
have the same q. We shall refer to the subspace spanned 
by these configurations as the projected space. For small 
9 the sin 9 term of the Hamiltonian mixes different config- 
urations within the projected space and also mixes states 
in the projected space with those outside. However, since 
there is an energy gap separating the projected space 
from the rest of the Hilbert space, the effective low-energy 
Hamiltonian, derived in degenerate perturbation theory 
to lowest order in sin#, is identical to the projection of 
H on the projective space. The Hamiltonian, with the 
ground state energy at 9 = subtracted, can be written 
as 

H = sin9V V. (6) 

i 

where V projects states into the projected space, Pij 
exchanges qi and qj, and the factor sin# defines simply 
an overall energy scale. 

The projected Hamiltonian is symmetric under the 
same discrete group as the full Hamiltonian. However it 
has two additional continuous U(l) symmetries. These 
symmetries are generated by the conserved charges N q , 
which measure the number of sites with the given value 
of q. Due to the constraint J2q=o = ^ omv two of the 
charges are independent. 
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The spectrum of Eq. ((6j) is invariant upon a global 
sign reversal of the Hamiltonian. The way to show this 
is to consider a division of the projected space into invari- 
ant subspaces, where each of these is spanned by tensor 
products of q-basis vectors, all having the same num- 
bers iVb,i,2- The states within such a subspace can all 
be mapped into each others by a permutation of the set 
of q values. A further subdivision is achieved by collect- 
ing all states that are connected by even permutations 
in one group. This defines two smaller subspaces, with 
basis vectors that are interconnected by odd permuta- 
tions. One can define a unitary operator T z that takes 
the value 1 in one of the subspaces and —1 in the other. 
Since the Hamiltonian ((6]) is the sum of transposition 
operators, it only has non-vanishing matrix elements be- 
tween states with different values for T z . This means that 
H anticommutcs with T z , hence the eigenspectrum of H 
is symmetric about E = 0. Upon 9 — > —9 the Hamilto- 
nian in Eq. ^ reverses sign. However the eigenspectrum 
remains unchanged. 

Now let us come back to the full Hamiltonian. Near 
9 = the low-lying energy spectrum changes linearly 
with sin#. Therefore, there is a massive crossing of en- 
ergy levels, with no level repulsion, at 9 = 0. In partic- 
ular the ground state changes abruptly, consistent with 
a first order transition at 9 = in the thermodynamic 
limit [L — > oo). 

The above discussion is in reality not restricted to m = 
3, but applies to all odd integer m. Thus, to lowest order 
in the deviation from 9 = 0, the low-energy Hamiltonian, 
for m = 2n+l, can be written as a projected Hamiltonian 
of the same form as Eq. ([5]) . Similarly one can show that 
there are now m — 1 conserved quantities, specified by 
the numbers No, ... , N m -x. The spectrum is inverted in 
the same way when 9 — > —9, and consequently, for odd 
integer m, there is at 9 = a first order phase transition. 
Due to the duality symmetry there is precisely the same 
type of phase transition at point 9 = tt/2. 



C. The critical region < 6 < n/2 

We now provide evidence for the two phases, one in 
the interval < 9 < tt/2 and the other in the interval 
— O.Itt < 9 < (tt/2 < 9 < 0.6 tt), which are connected 
by the first order phase transition at 9 = (tt/2), to be 
gaplcss. 

To determine the properties of the model in the region 
< 9 < tt/2, we have calculated the entanglement en- 
tropy by DMRG^. In Fig. @] we present a plot of the en- 
tanglement entropy Sa for a subsystem (subrcgion) A of 
length £ = L/2 for various different values of 9 in the in- 
terval < 9 < tt/2. In Fig.mja), 9 is fixed (tan0 = 1/5), 
while the number D of states kept during the DMRG 
calculation is varied. The plots show that for L = 400 
the value obtained for D = 1000 has almost converged. 
In Fig.rjjb), we compare the results for various angles 9, 
keeping D = 1000 states in each case. 
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FIG. 4: (a) The entanglement entropy associated with 
tan 9 = 1/5 as a function of InL under open boundary condi- 
tion. Different curves are distinguished by the different num- 
ber of kept states D in the DMRG calculation. The fitting 
curves are constructed using c = 2. (b) A comparison of the 
entanglement entropy for tan (9 = 1, 1/4, 1/5, 1/10 calculated 
with D = 1000. 



The scaling of the entanglement entropy with system 
size, shown in Fig. 01 is linear in InL, which indicates 
that for these values of 9 the model is described by 
a conformal field theory (CFT). For (l+l)-dimensional 
CFTs, with open boundary condition, the entanglement 
entropy Sa of a subsystem (subregion) A of length I, is 
S A = (c/6) log((2L/7ra) sin^/L)) + g + cq^L where 
c is the central charge associated with the critical be- 
haviour, a is a length scale, g is the boundary entropy, 
and Co a non- universal constant. In our calculation,we 
fix I to L/2, so that Sa — (c/6) log(L) for large L. 

From the numerical results, we conclude that for all 
values of 9 studied the system exhibits conformal invari- 
ance (hence is gaplcss) with a central charge c = 2. This 
value of the central charge can be understood by noting 
that for small 9 there are two conserved charges, each 
generating an emerging U(l) symmetry. For a similar 
gapless phase with c = 2, realized in interacting boson 
systems on the three-leg ladder at one-third filling, see 
Refs. [Mil 

Entanglement entropy plots for the region — 0.17T < 
9 < and 7r/2 < 9 < 0.67T show the same picture, and 
here too, we obtain c w 2 indicating gaplcssncss (see 
Fig. 5). The fact that the central charge is the same for 
these two 9 regimes can be understood as a result of the 
invariance of the spectrum of Eq. © upon reversing the 
sign of the projected Hamiltonian. 

To give further evidence that the system is critical, we 
also calculate the entanglement spectrum by means of 
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FIG. 5: Same types of plots as Fig. [4] here with values of i 
in the interval 7r/2 < < 0.67T. 



DMRG. In particular, we consider the eigenvalues {A;}of 
the density matrix. It was shown ir>i£ that the mean 
number of eigenvalues larger than a given A, denoted by 
n(A) can be calculated from CFT, with the result 



n(A) - / (2v / ^n(A max /A)) 



(7) 



where Jfe(x) is the modified Bessel function of the first 
kind, b = — ln(A max ) and A ma x is the largest eigenvalue 
of the reduced density matrix. 

In Fig. |5] we plot the distribution of the Aj for two 
points in the gapless region. Concretely, we plot the 
value of the ith eigenvalue, followingi£. The data agrees 
with the CFT result asymptotically, which confirms that 
the gapless phase can be described by CFT. In principle, 
it should be possible to extract more information about 
which CFT describes our system from the distribution of 
the small eigenvalues of the reduced density matrix, but 
this is a difficult talks, which we leave for future investi- 
gations. 

Numerical studies of energy gaps between low lying 
levels confirm the picture of gapless phases in the inter- 
vals referred to above. The DMRG evaluations of the 
gaps show a clear 1/L scaling behavior, consistent with a 
conformal field theory limit when L — > oo. As we already 
have pointed out, the numerical plots of the 9 dependence 
of the ground state energy show a discontinuity in the 
derivative at the points 8 = and 6 = n/2, consistent 
with a first order phase transition at these points. 

The space of c = 2 conformal field theories is rather 
extensive, see Ref. [TJI for an overview of possible theories. 
In the past, in order to identify the low energy conformal 
field theories of lattice models, it has been very instruc- 
tive to study the models under periodic boundary condi- 



FIG. 6: The distribution of the eigenvalues of the density 
matrix A. The length of the chain is 102 and the number of 
states kept in the DMRG is 1000. The blue solid line is the 
prediction of CFT. The step structure in the calculated data 
is due to the degeneracies of the eigenvalue of reduced density 
matrix. 



tions, and to calculate the energy spectra as functions of 
the various momenta. This can give crucial information 
in identifying the type of critical behavior. For the torus 




FIG. 7: The momentum resolved spectrum at 6 = 7r/4, 
for L = 18. The blue plusses correspond to the singlets 
(Q,Q) = (0,0); the red crosses to the quadruplets (Q,Q) = 
(1, 1); (1, 2); (2, 1); (2, 2) and the blue cross, red plus combina- 
tion to (Q,Q) = (0,1); (0,2); (1,0); (2,0). 

model we give in Fig. [7] the energy spectra, resolved with 
respect to momentum, for 6 = n/A, and for system sizes 
L = 12 and L = 18. 

From conformal field theory, it follows that if a system 
is critical, the energies of the low lying states are given 
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by 



27T7J C 

E = E L+— (h h + h K -—) 



(8) 



In this equation, Eo is a constant energy per site, v is a 
characteristic velocity, c the central charge, and h^, /ir 
the left and right scaling dimensions of the associated 
field in the CFT. In our case we know from the DMRG 
results that c = 2. Furthermore, from studies of the L 
dependence of the ground state energy, which has /il = 
hji = 0, the velocity parameter v has been determined. 
In FigJT] the energies are shifted such that the ground 
state has zero energy, and the energies are rescaled in 
units such that 2nv/L = 1. 

Scaling operators in CFT are grouped in terms of pri- 
mary fields, with total scaling dimension + = A p , 
and their descendants, with /il + hu = A p + n, where n 
is a positive integer. When the low-energy spectrum is 
presented in the form shown in Fig. [3 one should in prin- 
ciple, be able to identify the CFT through the position 
A p of the primary fields, with their associated towers of 
descendent state. However, in our case we have not been 
able to make such an identification. In part that is due 
to the wide range of possibilities for c = 2, as discussed 
in Ref. [l9|, but also due to the effects finite size effects for 
the systems we are able to perform exact diagonalization. 
Thus, the low energy spectra show several characteristics 
of a critical system, but we have not been able to con- 
clusively identify which c = 2 conformal field theory that 
describes the critical behavior. 



D. The phase transitions at 9 sa 0.67T and 6 
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The numerics indicates that a phase transition to a 
gapped phase, with a threefold degenerate ground state, 
takes place at 9 rj — 0.17T, with a similar transition at 9 ~ 
0.67T. These two gapped phases, which are connected by 
the duality transformation seem to cover the parameter 
range 0.6-7T < 9 < 1.9-7T. We will come back to these 
gapped phases in the next subsection. 

Fig. [Ha) shows the numerical estimate of the loca- 
tion of the critical point between the gapless phase and 
gapped phase in the vicinity of 9 m Q.6ir. In the plot the 
critical point is marked by the crossing between two exci- 
tation energies for different values of L. The curve "gap 
in the Q = sector" is the energy difference between the 
ground state and the lowest excited state in the Q = 
sector, while the curve >l E g (Q = I) — E g (Q = 0)" gives 
the gap between the lowest energy state in the Q — 1 sec- 
tor and the ground state. In the gapped phase the excited 
Q = state will merge with the ground state in the limit 
L — > oo, but in the gapless phase it lies above the lowest 
energy state in the Q = 1 sector. Thus a crossing be- 
tween the two excited states takes place at a point which 
moves towards the phase transition point as L — > oo. The 
plot shows a convergence of the crossing towards a point 
slightly below 9 — 0.617T. A plot of the second derivative 



o 

S-H 

-3 

>» 

S-H 

a 
w 



0.20 


—0—12 gap in Q=0 sector 


(a) 




1 


0.16 


. — 12Eg(Q=l)-Eg(Q=0) 

— - 18 ; ^ - 




—a 




0.12 










0.08 
0.04 




>- c 

3 ™ — . 




1 " 



0.600 



0.608 0.616 

9/71 



0.624 




0.615 0.620 0.625 0.630 0.635 0.640 
G/tc 

FIG. 8: (a) The crossing between the ground state in Q = 1 
and the first excited state in the Q = sector. Different color 
indicates different lattice size, (b) The second derivative of 
ground state energy with respect to 8. 



of the ground state energy through this point, as shown 
in Fig.tSJb), indicates that a continuous phase transition 
is taking place. It is worthy to note the continuous phase 
transitions in question link a gapless phase and a gapped 
phase. This is different from usual quantum phase tran- 
sitions where the phases on both sides are gapped. It is 
similar to the Kostcrlitz-Thouless phase transition except 
the central charge is different. 



E. The gapped phases 

We now focus our attention on the gapped phases of 
the m = 3 quantum torus model. A first gapped phase 
extends from the phase transition near 9 w 0.67T to 9 = 
5ir/4. The second gapped phase is dual to the first one, 
and ranges from 9 = bn/A to 9 pa I.97T. We will discuss 
the phase transition between these two gapped phases in 
the following subsection. 

To understand the symmetry breaking in the gapped 
phase it is useful to consider the special points 9 = tt 
and 9 = 3ir/2. For these values of 9 the Hamiltonian 
becomes classical (i.e., there is no non-commuting oper- 
ators), and where the ground state is exactly threefold 
degenerate even for finite L. Away from these points 
there is finite-size lifting of the ground state degeneracy. 
However the latter decreases exponentially with increas- 
ing L. In Fig. EJa) the presence of a gap is demonstrated 
via the saturation of the entanglement entropy as a func- 
tion of InL. The exponential decay of the finite-size gap 
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between the states which evolve into the ground state is 
shown in Fig. [3(b). 

At the special points 9 = w the threefold degenerate 
ground state is spanned by the product states Yli=i W)> 
m = 1,2,3. As a consequence the correlation function 
(If® If?) is simply a constant, while the corresponding 
U v correlation function vanishes for i j. Even if this 
is a very special situation, a similar behavior of the cor- 
relations functions is seen in the full gapped phase, up 
to 9 = 57r/4. Thus, the U x correlation function is long 
range, while the U y correlation function decays expo- 
nentially with the distance between the two points i and 
j. This is shown for a particular value 9 = 1.2657T in 
Fig. [TU] for a system of length L = 102. At the point 
8 = 57r/4 there is an interchange between the correlations 
of U x and U v , as follows from the duality symmetry. The 
curve corresponding to 9 = 57r/4 in Fig. 1101 shows iden- 
tical long range correlations for U x and U v . This can be 
understood as due to the fact that the ground state for 
any finite L is a 50/50 superposition of two states, with 
correlations that are symmetric under the interchange 
x o y. 



F. The phase transition at 9 = 57r/4 

As a final point, we have examined the phase transi- 
tion at the point 9 = 5ir/4, at which the two gapped 
phases meet. In the vicinity of 9 = 57r/4 there are six en- 
ergy eigenstates that are important: two (Q, Q) = (0, 0) 



singlets, a (Q,Q) = (0,1); (0,2) doublet and another 
(Q,Q) = (1,0); (2,0) doublet. All the numerical evi- 
dences are consistent with the following picture in the 
thermodynamic limit. For 9 < 57r/4 the ground state 
is triply degenerate. The triplet involves one of the 
(Q,Q) = (0,0) singlet and the (Q,Q) = (0,1); (0,2) 
doublet. For 9 > 57r/4 the ground state is also triply 
degenerate. This time the triplet involves the other 
(Q,Q) = (0,0) singlet and the (Q,Q) = (1,0); (2,0) 
doublet. At 9 = 5ir/4 the two triplets cross, result- 
ing in a six fold degenerate ground state. For finite L 
the doublet (Q,Q) = (0,1); (0,2) precisely degenerate 
with (Q,Q) = (1,0); (2,0). This is because they form 
the four dimensional irreducible representation of the en- 
larged symmetry group at the selfdual point. In contrast 
the two (Q, Q) = (0, 0) singlets are slightly split due to 
avoided crossing caused by the finite system size. The 
above picture suggests a first order phase transition at 
the selfdual point in the thermodynamic limit. This is 
caused by the crossing of energy levels. The two phases 
are distinguished by different long range correlations for 
either U x or U v , as previously discussed. 

The character of the phase transition is further illus- 
trated by the plots in Figs. El Hand EI] In Fig. El we 
display the behavior of the gap at the transition point 9 = 
57r/4. The ground state is a (Q,Q) = (0,0) singlet, and 
Ai denotes the energy difference with the first excited 
'level', which is the (Q,Q) = (0, 1); (0, 2); (1, 0); (1, 0) 
quadruplet (see the discussion above). A 2 denotes the 
energy difference between the ground state and the sec- 
ond excited 'level', again a (Q,Q) = (0,0) singlet. The 
energy difference Ai and A 2 are given in Fig. [TT] for vari- 
ous values of D, the number of states kept in the DMRG, 
clearly showing that the results have converged as a func- 
tion of D. Both the energies of the first and second ex- 
cited levels decay exponentially with system size, which 
implies a six-fold degenerate ground state, because other 
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FIG. 10: Logarithmic plot of the U x and U y correlations at 
= 57t/4 and 6 = 1.2657T. For 6 = 5tv/4 the two curves fall 
on the top of each other, d is the distance between the two 
points for which the correlation is calculated. The system size 
used in the calculation is L — 102. 
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level (a quadruplet, see the main text) and the ground state 
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levels have a finite gap to these six degenerate states in 
the thermodynamic limit, see also Fig. [T2] 

Panel (a) of Fig. [12] shows the behavior of the energy 
difference between the two lowest lying levels in the Q = 
1 sector. The latter does not merge with the degenerate 
ground state in the thermodynamic limit. The plot shows 
that a small energy gap remains as L — > 00. 

In panel (b), the first derivative of the ground state 
energy is plotted function of 9. It shows a sharp 
drop at 9 = 57r/4, which is expected to become infinitely 
sharp in the limit L — ¥ 00, because transition happens 
at an isolated point between two gapped phases. The 
sharp change in the derivative of the ground state energy 
is consistent with the picture of the phase transition as 
caused by the crossing of the two lowest energy levels. 

In Fig. [131 we show the results of a detailed study 
of the behavior of the entanglement entropy, near the 
point 9 = 57r/4. The plot shows that for system size up 
to about L = 256, the entanglement entropy behaves as 
expected for a critical point, and is consistent with c = 2. 
For larger system sizes, the entanglement entropy flattens 
of, and finally crosses over to the (lower) value obtained 
via the MPS method. This last crossing over happens 
at larger system size, if the number of states kept in the 
DMRG calculation is increased, although for D = 2000, 
the numerics seems to have converged. 

The drop in the entanglement entropy for large L, 
at 9 ~ 57r/4, can be understood in terms of the level- 
crossing picture as follows: There are two states with 
{QiQ) = (0,0) participating in the crossing; let us de- 
note them \ip-) and |V>+), where \tp—) (|"0+)) is a ground 
state for 9 < 5tt/4 (9 > 57r/4). They are not orthogonal 
but the overlap is exponentially suppressed for large L. 
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FIG. 13: The entanglement entropy, close to the point 
57t/4, as obtained from DMRG and the MPS method. 



For finite L, the true ground state is a mixture of |^/>-} 
and 1^+), which at 9 = 57r/4 is an equal-superposition 
thereof. Away from this point the ground state rapidly 
rotates into either or 1^+), depending on whether 
9 is reduced or increased from 9 = 5ir/4. This rotation 
is more rapid the larger L is so that for L — > 00 the 
crossing between \ip-} and IV'-i-) becomes sharp. Assume 
we choose 9 slightly smaller than 57t/4. For small L the 
ground state is essentially the 50 percent mixed state, 
but when L increases at some point it rapidly changes 



10 



to The entanglement entropy of each of the two 

states and \ip-) clearly is smaller than that of the 
superposition of the two. Due to the symmetry between 
and at 9 = 5tt/A, and the rapid decay of the 
overlap between them as L — > oo, we can in fact estimate 
the drop in entanglement entropy. Since each of the two 
states in the superposition gives equal contributions to 
the entanglement entropy the drop should be close to 
In 2. The numerical value found for the curves plotted in 
Fig. [T3] is in fact very close to In 2. 

We have conclusively shown that the gap at the tran- 
sition point at 9 = 5tt /4 does not close, and that that 
system exhibits long-range order. It is intriguing to ob- 
serve that the scaling of the entanglement entropy shows 
behavior consistent with c = 2 critical behavior) up to 
fairly large system sizes of at least L = 200 (due to the 
rather large correlation length). Combined with exact di- 
agonalization results (which we performed up to L = 18, 
not shown), one could incorrectly be led to believe that 
the system is critical. It would be interesting to investi- 
gate if the system can be (fine-) tuned to become criti- 
cal, by perturbing away from the gapped 9 = 57r/4 point. 
One could for instance allow complex amplitudes in the 
Hamiltonian (|5|), or mixed terms such as UfU^ +1 . 

V. CONCLUSIONS 

To summarize, we have presented a class of one- 
dimensional lattice models, called the quantum torus 
chain. These models depend on an integer-valued pa- 
rameter m (the number of magnetic flux quanta piercing 
through each torus) which controls the presence/absence 
of frustration. In addition there is a continuous parame- 
ter 9 which controls the ratio of the non-commuting terms 
in the Hamiltonian. The models show a variety of gap- 
less and gapped phases. In all cases the gapped phases 
break symmetries of the Hamiltonian. We characterize 
these gapped phases based on the properties of specific 
points in parameter space, where the Hamiltonian is ex- 
actly solvable. For the gapless phases, we have used 
DMRG and matrix product state methods to evaluate 
the entanglement entropy and hence deduce the central 
charge. 

Based on these analysis we have concluded that while 
for even m the system is gcncrically gapped with iso- 
lated quantum critical points separating different gapped 
phases. For odd m there exists extended regions in the 
parameter 9 where the system is gapless. The above gen- 
eral conclusion is motivated by the mapping of the m = 2 



model on the spin 1/2 XZ model, and the detailed nu- 
merical studies of the m = 3 model. For the latter the 
quantum phase transitions between two gapless or two 
gapped phases are first order, while the transitions be- 
tween gapless and gapped phase are continuous. Consis- 
tent with these conclusions we have numerical evidence 
from studies of the m — A model that this is also gapped 
in the first quadrant < 9 < tt/2. 

Our study raises a series of open questions. At the 
moment we have not been able to pin down the conformal 
field theory for the c = 2 gapless phases of the m = 3 
model. Of course gapless phases have precedents in one 
dimension (e.g. the Luttinger liquid^). However gapless 
phases in more than one dimension are very unusual. Our 
model of course can be defined in any spatial dimension. 

Without going into much detail, we can say that for 
m even, there is no frustration and at the special angles 
9 = 0, 7r/2, 7r, 37r/2, the ground states are direct product 
states, with a finite degeneracy. Thus at least near these 
points, we expect the system to be gapped. Whether or 
not these gapped phases extend over the full parameter 
space, leaving only critical points as is the case in our 
1-d model needs further investigation. For odd m, there 
is also quantum frustration in higher dimensions. It is 
interesting to ask whether the gapless phase still exists 
in higher dimensions. 

There are several ways to extend our model. One could 
try to include terms which mimic a magnetic field in or- 
dinary spin chains. In our model, the closest analog of a 
magnetic field would be a one-site term, which breaks the 
symmetry of the model, by favoring (say) one of the pos- 
sible states. More phase transitions should be expected 
in this enlarged parameter space, because a strong sym- 
metry breaking term will simply 'polarize the degrees of 
freedom', leading to ordered phases. 

The general properties we find in the this work are 
consistent with the results discussed by Chen et alX for 
systems where the internal symmetry group is repre- 
sented locally on each site by a projective representation. 
It will be interesting to consider generalizations of the 
torus model where the symmetry is represented linearly 
locally. In this case symmetry protected topological 
states become possible. Finally it is interesting to ask 
whether it is possible to describe the gapless phase in 
terms of a discrete non-linear a model with a topological 
term. We will leave these and other related questions as 
interesting projects for future investigation. 
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